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^ ■ Abstract 

In the large U limit, the ground state of the half-filled, nearest-neighbor Hub- 
^ . bard model on the triangular lattice is the three-sublattice antiferromagnet. 

. In sharp contrast with the square-lattice case, where transverse spin-waves 

. and charge excitations remain decoupled to all orders in t/U , it is shown that 

0^ ■ beyond leading order in t/U the three Goldstone modes on the triangular lat- 



tice are a linear combination of spin and charge. This leads to non-vanishing 



S 

. conductivity at any finite frequency, even though the magnet remains insu- 

O ■ fating at zero frequency. More generally, non-collinear spin order should lead 



to such gapless insulating behavior. 
PACS numbers: 71.30.+h,75.30.Ds, 72.20.-i,75.30.Fv 

In band theory, insulating behavior is due to the existence of an energy gap in the single- 
particle excitation spectrum. It is well known, however, that interactions neglected in band 
theory may lead to insulating properties, a point first discussed by Mott and, in 1964, by 
Kohn who showed that the existence of an energy gap is sufficient but not necessary to 
have an insulator An insulator is thus best defined by a vanishing zero-temperature DC 
conductivity. In the usual examples of Mott insulators, however, there is generally a gap 
not only in the single-particle excitations, but also in the frequency dependent conductivity 
which is related to excited states at constant number of particles (particle-hole excited states 
instead of single-particle excitations). In this letter, we show that when interactions lead to 



magnetic order in itinerant spin 1/2 systems, the resulting Mott insulator generally has no 
gap in the frequency dependent conductivity: This is the behavior of a gapless insulator. 
It is only when the magnetic order is collinear antiferromagnetism on hypercubic lattices 
that one can demonstrate the existence of a gap in the conductivity. In the general case 
that we consider in detail, the mixed spin and charge character [0] of the finite q collective 
(Goldstone) modes leads to the vanishing of the gap in the conductivity for a finite range of 
interaction strengths, despite the fact that the single-particle excitations do have a gap. 
We start from the general one-band Hubbard model 

H = - ^ij {4aCja + c]„Ci^) + ^ Ui^n^ , (l) 

where the first sum is over pairs of sites, tij is the hopping integral and U the on-site Coulomb 
interaction. We then define 

^"(0 = ^Ec.X,c./3, (2) 

a dimensionless four-component spin and charge operator. (To be explicit, we let the matrix 
index /i take the values p,x,y,z, and we define a'^ as the unit matrix.) We consider the 
general class of magnetically ordered states with planar spiral order 

{Sn + ^{Sn = Se^"^-^^ , (3) 

where S is the order parameter representing the average moment on each site. Without loss 
of generality, we take the spiral in the x — z plane. Well known examples of such states 
include the collinear antiferromagnet on the square lattice, Q = 7r/ax + n/a z, and the 
three-sublattice antiferromagnet occurring on the frustrated triangular lattice, a 120*^ spiral 
with Q = A7c/3a X. As always when working with planar spiral order, the analysis is greatly 
simplified by using a rotating (orthonormal) basis in which the quantization axis at every site 
points in the same direction as the average spin density of Eq.(^. In this rotating frame, the 
single-particle Green's function Gafsii^j] t) = — (Tcja(T)cj^(O)) becomes, in matrix notation, 

G{iJ; iun) = T} G{i,j; iujn) Tj , (4) 
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where the rotation matrix in spin space Ti = e-^(Q-^')^''/2 (with ay the usual Pauh matrix) 
depends on the site i. In this reference frame, the Green's function matrix G has the full 
underlying lattice periodicity, even for incommensurate spiral order. In the Hartree-Fock 
approximation, we find 

^' ia;„ + At-E+(k) iu;„ + /x-£;_(k) ' ^' 

where 



^ E{k) E(k) 

V +^E(k) E(k) y 



with the definitions 



eo(k) = 5^ cos k • (R. - R^) 

for the single-particle dispersion in the paramagnetic phase ( is the number of sites) , and 

e± = eo(k±Q/2) r^{\,) ^ ^-±^1^ 
e(k) = + 1 ^ = US ■ (6) 

E{\i) = <^?72(k) + A2 £;±(k) = e(k) ± £;(k) 
The gap equation is given by the self-consistency requirement, 

u ^ [/(£;_(k))-/(E+(k))] _ 

2Ar^ £;(k) ' ^ ' 

where f{x) is the Fermi function and the wave vector k spans the entire paramagnetic 
Brillouin zone of the crystal. The chemical potential is determined by number conservation, 

^E[/(^-(k)) + /(i^+(k))] = l. 
k 

Proceeding to the collective excitations in the ordered state, we define the matrix response 
function 

x'^'ihr.T) ^ - {TS^{i-T)s^{j-Q)) + {s>^{i)) {S'^ij)) . 

X is obtained in the Generalized Random-Phase Approximation (GRPA) by the usual sum- 
mation of bubble and ladder diagrams. Since the GRPA is conservative, going beyond this 



approximation should not change the quahtative aspects of ours results. In the rotating 
frame where x^^ihJ'i ~^ X^^i'i ~ J] t) , the matrix GRPA equation takes the form 

X(q, = X °(q, i^n) + 2Ux °(q, ^f^„)rx(q, i^n) , (8) 

where F is a diagonal matrix with F^^ = 1, F^^ = F^^ = F^^ = — 1. 

At T = OK, which we will consider from now on, the retarded zeroth-order matrix 
susceptibility x "(q, + i6) is given by, 

/ 



^)A- ^A(ii;)A_ (i^ + ^)A+ A(^ + i.)A 
^A{^)A_ {^±^)A- .A(^ + i.)A, (| + ^) A 
+ ^ ) -.A (1 + i) A, (A!±||;-^) A_ -A (^) A_ 

A(^ + i.)A, .(1 + ^)A^ -A(^)A_ (-A!^^). 

(9) 

where unprimed functions are to be evaluated at k+q/2 and the primed functions at k— q/2, 
and where A± = 7_ ± 7_(_ with 

" u; + iS~{e-e')±{E + E') ' 
The response function in the rotating frame x can be related to the laboratory response 
X but the expressions are lengthy [Q. For the density response however, the relation simply 
is 

^'"'(q, q', uJ + i5)= r^lq, uj + i5)5^,^ . (11) 

The other response functions involve linear combinations of x at different wave-vectors. The 
poles of Xi which of course coincide with those of Xi give the position of the collective modes. 

For hypercubic lattices in arbitrary dimension larger than one, symmetry arguments can 
be used to show that whatever the hopping matrix tij and the value of f/, as long as there is 
long-range collinear antiferromagnetic order, the response matrix x is block diagonal. There 
are then two Goldstone modes with purely transverse x — y character. For the nearest- 
neighbor model at half-filling, it is well known 0] that in the limit t <^ U, the above 



approach correctly reproduces the standard spin wave result of the Heisenberg model with 
J ~ At^/U. For next-nearest neighbor hopping on the square lattice, the resulting spin 
Hamiltonian in the large U limit is frustrated. Nevertheless, it follows from above that as 
long as there is long-range antiferromagnetic order, there are two Goldstone modes with 
purely transverse x — y character, as in the simpler non-frustrated case. 

From now on, we restrict ourselves to the nearest-neighbor triangular lattice at half- 
filling. On this lattice, the frustrated antiferromagnetism can give rise to non-collinear 
magnetic order, a subject of current interest in itself Much work has been done to obtain 
the magnetic phase diagram as the electrons are allowed to be more and more itinerant {i.e. 
as the ratio U/t decreases) [||. Krishnamurthy et al. and Jayaprakash et al. P have 
shown that, at half- filling, as U is increased the Hartree-Fock ground-state evolves from a 
paramagnetic metal to a metallic incommensurate spiral Spin-Density- Wave (SDW) then to 
a commensurate linear SDW with indirect gap and finally, at higher values of f/, into a spiral 
SDW insulator with a three-sublattice 120° twist between spins on neighboring sites. This 
spiral (or equivalently helical) SDW state is the well known ground state of the Heisenberg 
model (with J ~ At^/U). We call this state the 120° spiral SDW phase. 

We go beyond previous studies by accounting for the collective excitations of the 
120° spiral SDW phase. The electronic dispersion in the absence of interaction for the 
nearest-neighbor model on the triangular lattice is given by eo(k) = — 2t(cos(A;^.a) + 
2 cos(A;^a/2) cos(v^fcyCt/2)). In the paramagnetic Brillouin zone, the single band eo(k) is 
split into the two subbands -E±(k) by the presence of spiral order, as can be seen from 
Eqs.(|^),(^. The direct single-particle energy gap is 2 A but the indirect gap Ejg = 
Min[£'+(k)] — Max[£'_(k)] can be lower than this value. When the indirect gap vanishes, a 
transition to a metallic state occurs. This transition takes place at a value of t/U ~ 0.195, 
above which the order parameter S decreases very rapidly. However, the collective modes 
destabilize the insulating 120° spiral SDW phase at t/U ~ 0.146, well before the above metal- 
lic transition is reached. In the following, we restrict ourselves to the regime at half-filling 
where the indirect gap is positive and only the lowest energy band is filled. 
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Expanding x°(q, tu) Eq.(|^) to second order in uj/U,t/U, and using the gap equation to 
eliminate 5* in favor of f/, we find after tedious algebra that to leading order in t/U , the poles 
co'(q) of the response function Eq.(|) correctly reproduce the standard spin wave result of 
the Heisenberg model with J = At^ /U [jl0|. If q is restricted to the magnetic Brillouin zone, 
the three spin-wave branches (Goldstone modes) are given by co'(q), ^^(q + Q),c<j(q— Q). 

As the ratio t/U increases and the electrons become more itinerant, higher-order terms 
mt/U cannot be neglected and the coupling between charge and the three spin components 
becomes important. The matrix x(q, cu) is no-longer block diagonal and all response func- 
tions then share the same poles, although with different weights. Three of the poles are still 
Goldstone modes with a vanishing frequency at q = in the magnetic Brillouin zone. In 
Fig. 1 we plot the imaginary part of the response functions x'''' ^ind fo^' various values 
of q. We plot in Fig. 2 the dispersion relation of the collective modes obtained by following 
the peak of the response function as q is varied along high-symmetry directions. 

It is clear from Fig. 1 that the coupling between spin and charge is very small so that 
the spin waves are only very slightly modified by the charge fluctuations. The intensity of 
this coupling increases with t/U but there is only a narrow range of t/U over which the 
120° spiral SDW is stable. Indeed, at t/U ~ 0.146 one of the collective modes softens, 
destabilizing the 120° spiral SDW ground state at a smaller value than the t/U = 0.16 
found in Ref. |^ using only self-consistent Hartree-Fock solutions. 

Because the collective modes with coupled spin-wave and charge extend to zero-frequency, 
it is natural to ask whether the system could becomes conducting despite the existence 
of a gap in the single-particle excitations. The appropriate definition of an insulator is 
that it has a vanishing Drude weight. |l|,0 This statement is equivalent to saying that 



the coefficient of the delta function S{uj) in the zero-temperature DC conductivity anc = 
lim^^o liKifj^o '^(^) q) vanishes. We obtained this quantity from the appropriate limit of 

aiu, q) = ^e'^-^^^±^x''{u; + q, q) . (12) 
In this expression, one has to take the part of x''''{'^ + ^"^j q) =X^^{'^ + q) which is 



irreducible with respect to the interaction. This ehminates the effect of screening. Ex- 
panding to second order in q, cj, we found after lengthy algebra that the charge-response 
function decreases faster than for g ^ so that the Goldstone-mode contributing to 
o"(ci;, q) has effectively zero weight at q ^ 0. The same conclusion is reached by calculat- 
ing the transverse current- current response function in the GRPA. It then follows that the 
system remains insulating even though there is no gap in the charged collective excitations. 
This is not so surprising since at q = the Goldstone modes restore rotational invariance: 
hence they have the same parity as the ground state so that matrix elements of the current 
operator between the q = modes and the ground state vanish. 

At finite wave vector and frequency, however, the conductivity is finite and, as shown in 
Fig. 3, the absorption described by the real part is exactly at the Goldstone-mode position, 
as expected. We conjecture that the insulator can be changed into a bad conductor by 
impurities. Indeed, impurities transform the single delta function 5{uj) appearing in the DC 
conductivity of perfect metals into a Lorentzian-like response. In our case, the breaking 
of crystal momentum conservation by impurities would probably, by a similar mechanism, 
broaden the finite q resonances to make the g = conductivity finite everywhere in the 
single-particle gap, including zero frequency. The question is difficult because of possible 
Anderson localization. 

In conclusion, the triangular three-sublattice antiferromagnet described by the large U 
half-filled one-band Hubbard model provides an example of a gapless Mott insulator. In 
this system, the Goldstone modes have mixed spin and charge character, leading to the 
disappearance of the gap in the conductivity, despite a vanishing DC conductivity and the 
existence of a gap in the single-particle excitations. This phenomenon is generic for itinerant 
magnets with spiral order. An exception is when the ordering wave vector corresponds to 
coUinear antiferromagnets. 
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FIGURES 

FIG. 1. Imaginary part of the density and (inset) spin response functions at 

t/U = 0.1. From left to right the peaks, broadened by a smah imaginary part, correspond to 
qy = 0, = 0.1, 0.15, 0.25 (in units of 27r/a) 

FIG. 2. Dispersion relation of the mixed spin and charge collective modes along the edges of 
the irreducible Brillouin zone of the triangular lattice. The full, dashed and dot-dashed curves are 
respectively for t/U = 0.1,0.12,0.146. 

FIG. 3. Real part of the conductivity defined in Eq. (13) at t/U = 0.1 and for the wave vectors 
Qy = 0; = O-lj 0.15, 0.25 (in units of 27r/a) represented by the full, dashed and dot-dashed curves 
respectively. 



10 



0.08 




wave vector 




0.00 0.02 0.04 0.06 

co/U 



